Quasi one-dimensional magnetism driven by unusual orbital ordering in CuSb2 06 
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Essentially all undoped cuprates exhibit a quasi-planar, fourfold Cu-0 coordination responsible 
for the magnetically active antibonding 3d^2_y2-\ike state. Here, we present an electronic structure 
study for CuSb206 that reveals, in contrast, a half-hlled 3d^^2_^2 orbital. This hitherto unobserved 
ground state originates from a competition of in- and out-of-plaquette orbitals where the strong 
Coulomb repulsion drives the surprising and unique orbital ordering. This, in turn, gives rise to an 
unexpected quasi one-dimensional magnetic behavior. Our results provide a consistent explanation 
of recent thermodynamical and neutron diffraction measurements. 
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Low dimensional systems have always been of funda- 
mental interest to both experimentalists and theorists for 
their peculiar magnetic properties. In recent years, espe- 
cially boosted by the discovery of the high- Tc supercon- 
ductivity in cuprates, low dimensional magnets related 
to this family of compounds have been widely studied. 
Low dimensional cuprates exhibit a large variety of ex- 
otic physical properties. This variety is a result of the 
complex interplay of different interactions; mainly cova- 
lency, ligand-fields and strong correlation effects. 

A nearly universal component of cuprate systems is a 
strongly elongated CuOg-octahedronllJ wherein the ex- 
otic behavior finds its origin in the deceivingly simple 
planar, half-filled CU-O4 orbital lying in its basal plane. 
Here, Cu 3d and O 2p states form a well separated, half 
filled pd ~ a molecular-like orbital of x'^ — symme- 
try that is magnetically active. In general, the low-lying 
magnetic excitations, especially the magnetic ground 
state can be well described within this subspace of the 
whole Hilbert space. This way, the low-dimensional mag- 
netic properties of a compound depend crucially on the 
arrangements of these CU-O4 units forming different net- 
works, for instance (i) quasi ID chains from isolated_(i.e. 
Sr2Cu(P04)2@, i) edge-sharing (i.e. LiaCuOai, l5|)or 
corner-sharing (i.e. Sr2Cu03[l,[3) CUO4 plaquettes, (ii) 
quasi 2D square lattices in La2Cu04[l|, etc. 

In contrast to the vast majority of cuprates with a 
quasi-planar Cu(II) four-fold coordination, copper dianti- 
monate CuSb206 - first synthesized in early 1940's[9i 
- exhibits a nearly octahedral local environment of the 
Cu^"*" ions. The compound is green in color [§, in- 
dicating insulating behavior with a gap size typical for 
cuprates. Although the compound undergoes a second 
order phase transition below 380 K from a tetragonal 
to a monoclinically distorted structure, [ll[ this quasi- 
octahedral local environment is basically unchanged. 
The magnetic cation snblattice is that of the K2NiF4 
structure type, which includes many classic examples of 
square lattices exhibiting 2D antiferromagnetism ( AFM) , 
for instance the above mentioned La2Cu04. Surprisingly, 
susceptibility measurements done on both powder and 
single crystal samples of CuSb206 fit extremly well over 



a large temperature range to a nearest-neighbor-only S 
~ i Heisenberg ID model with an exchan ge c onstant 
ranging from -86 K to -98 K.[l3, [ll [13, \I£M, M, [13 
Furthermore, all low temperature susceptibility measure- 
ments show a sharp drop at 8.5K(ToI [13. [isl. [l^ [isl. [l6l . 
due to the onset of long range AFM ordering, [ij, [l^, [iq 
Thus, the contrast of the 2D lattice and the ID magnetic 
behavior in CuSb206 has been puzzling up to now. 

In this paper we present the combined results of the to- 
tal energy and tight-binding model (TBM) calculations, 
which indicate an hitherto unobserved ground state orig- 
inating from a competition between the in- and out-of- 
plaquette orbitals. Strong Coulomb correlation drives a 
surprising and unique orbital order in CuSb2 06, thereby 
leading to the strongly one-dimensional magnetic behav- 
ior with exchange integrals in good agreement with those 
deduced from experiments. Our results are consistent 
with recent neutron diffraction measurements. [l6j 

The tetragonal structure[ll| (a = 4.629 A , c = 9.288 
A) of CuSb2 06 with two formula units per unit cell is elu- 
cidated in Fig. [T] One interesting aspect of this system is 
the bond lengths, the out-of-plaquette Cu-0 bond length 
is shorter by 0.04 A compared to the in-plaquette Cu-0 
bond lengths. The CuOg octahedra are well separated 
from each other, with no edge-sharing or corner-sharing 
between adjacent CuOg octahedra. An angle (3 « 91.2° 
is the main distinguishing factor of the monoclinic (T < 
380K ) phase of CuSb206. Nevertheless, the calculated 
density of states (DOS) and band structure for both these 
phases are nearly identical. Therefore, all band structure 
calculations shown here have been performed using the 
tetragonal tri-rutile structure instead of the monoclinic 
one for the sake of simplicity. Using the more regular 
tetragonal structure is quite advantageous in calculating 
and understanding the tight-binding hopping parameters 
and in the construction of larger supercells to model spe- 
cific magnetic orderings. 

We performed band structure calculations within the 
local (spin) density approximation (L(S)DA) using the 
full-potential local orbital minimum basis code FFLO 
(version 5.00-18). [H [H In the scalar relativistic calcu- 
lations the exchange correlation potential of Perdew and 
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FIG. 1: (Color Online) The tetragonal trirutile structure of 
CuSb2 06. There are two Cu atoms per formula unit (two 
sublattices) . The CuOe octahedra in these 2 sublattices are 
rotated by 90 with respect to each other. Both the global 
{x,y,z) and the local (i, y, z) coordinate systems used in our 
calculations are depicted. Within the local coordinate sys- 
tem, the in-plaquette oxygens are highlighted in grey while 
the out-of-plaquette oxygens are in green. The Sb atoms are 
highligted in yellow. The in-plaquette Cu-O bondlength is 
0.04 Alarger than the out-of-plaquette Cu-O bondlength. 



Wang[21i] was used. Additionally, the strong Coulomb re- 
pulsion in the Cu id orbitals was taken into account in a 
mean field approximation by LSDA-I-C/ calculations [20j. 
using the around-mean-field double-counting scheme. 
Variations of Ud {J = leV) within the relevant physi- 
cal range {Ud = 5 to 9 eV) did not change the reported 
results significantly. 

The calculated non-magnetic LDA band structure, dis- 
played in Fig. [21 shows four bands crossing the Fermi level 
Ei?, i.e. two pairs of narrow and broad bands. Consider- 
ing the two formula units per cell, this is a very unusual 
result. In the standard cuprate scenario, only one well 
separated half-filled antibonding band per Cu^"*" is ex- 
pected. Projecting the orbital character onto the band 
structure in a local coordinate system wherein the local 
z-axis is perpendicular to the Cu-O plaquette-plain (see 
Fig. [T]) and points towards the apical oxygens, we find 
that the 'narrow' bands originate from the in-plaquette 
Cu dx2-y2 orbitals, while the 'broad' bands come from 
the out-of-plaquette Cu d^^2_f:2 orbital. This is differ- 
ent from the conventional cuprate scenario where only 
one band per Cu atom crosses the Fermi level and is 
comprised predominantly of the in-plaquette Cu 3d or- 
bital character. This unique feature stems from the fact 
that the CuOg octahedra are only slightly distorted in 
CuSb206, so that the cubic degeneracy for the Cg ligand- 
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FIG. 2: (Color Online) LDA band structure of the tetrag- 
onal CuSb206 along high-symmetry directions. All four Eg 
bands are partially filled. The 'narrow' band (open circle, 
blue) is entirely from the local Cu d^2_y2 orbital, while the 
'broad' band (filled circle, red) comes from the Cu d^^2_f2 
orbital. The inset shows the leading hopping terms of the 
tight-binding model (see text). 



field states are only slightly lifted. The energy difference 
between the c?j2_j^2 and d3j2_,~2 related band centers is 
about 0.3 eV only, to be compared with about 2 cV for 
'standard' cuprates. 

Because of the strong competition between the in- and 
out-of-plaquette orbitals, the effects of strong correla- 
tions (LSDA-I-C/) on the magnetic model are not imme- 
diately obvious. Let us first review the possible ways in 
which ID order can be invoked in this system. Scenario 
1: The interaction of the electrons is mediated via the 
Cu04 plaquette, ID chains will be formed along the [110] 
(and along [110] in the other sublattice). Scenario 2: 
The interaction of the electrons is mediated via the out- 
of-plaquette oxygen. This would in turn create ID chains 
along the [110] (along [110] in the other sublattice) direc- 
tion. These two competing scenarios in principle could 
lead to the same spin ordering pattern of ID antiferro- 
magnetic chains. Therefore strong correlations may drive 
the system to either one of these two scenarios. 

Firstly, to compare the microscopic magnetic interac- 
tions in CuSb206, we did a two-site one-band TBM fit, 
separately for the dj2_j^2 and d3j2_^2 band complexes. 
The respective 2x2 tight-binding matrix is given by. 



a + (3 S 
6 a + (3' 



with 

a = 
P - 

a' - 

5 - 



eo + 2h [Ci + Cy] + 2UCi + 4t5 [d + Cy] d 

2^2 [CxCy ~ SxSy] -\- 2t^ [CxCy -\- SxSy\ 
2^3 [CxCy — SxSy] + 2t2 [CxCy -\- SxSy] 

8tRCi;fCi,~,Ciz + StiCiz Cs sCij-. — CiiiCa,- 

2 2^ 2 2 2 2^ 2 2^ 
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where Cx = cos(fcja), Ci ^ = cos(fci|), Cs^ = cos(fc£^), 
Si = svD.{kxa) and the corresponding cyclic permuta- 
tions. The values of the transfer integrals obtained from 
the fit are collected in Table. HI with the hopping paths 
shown in the inset of Fig. [2l The hopping parameters 
obtained from the fit to the narrow dj2_jj2 band com- 
plex are all quite similar in magnitude, implying almost 
equal interaction strength along the three main hopping 
paths ti, t3, and t4 (see Table. HI This is indicative of a 
3D magnetic model and the energy scale of the transfer 
integrals are too small. This 3D model is incompati- 
ble with the experimental data.[l^ The fit to the broad 
d3f2_j:2 band produces a transfer integral (t2) which is 
very large compared to all others (see Tabic. This 
suggests a strong ID AFM interaction along the [110] 
direction, implying that the superexchange path Cu-0- 
0-Cu is along the [110] basal diagonal and is via the 
apical oxygens and NOT via the in-plaquette oxygens. 
This is in line with the band characters of the projec- 
tion mentioned above. This leads to almost perfect ID 
chains of antiferromagnetic Cu atoms along the [110] di- 
rection. Another interesting fact to notice is the differ- 
ence in the next-nearest- neighbor interaction along [110] 
{t2) and [110] {ts). The Cu-0 bondlength along [110] is 
shorter (2 A ) with a 180° bond angle as compared to the 
Cu-0 bondlength of 3.5 A along [110] and a 160° Cu-0- 
Cu bond angle. This structural feature also helps in the 
one dimensionality along [110]. The individual exchange 
constants calculated using J^- = 4t|^/Ue//, with a Ue// 
= 4.5 eVfil] are also collected in Table. D The sublat- 
tice coupling, connecting the two Cu sites in the unit 
cell is rather small, but not negligible, and so are the 
interchain couplings. This is very much consistent with 
the ID nature of the system. Though the plaquettes are 
quite isolated in CuSb2 06, they can interact in many dif- 
ferent ways and this in turn determines the magnitude of 
the exchange constants. But the unique geometry in this 
system gives a very large ID exchange (J^^*^ = 400K) 
as compared to other 2D cuprates. This value of J2 is 
too large (by a factor of four) due to possible ferromag- 
netic contributions like in other cuprate compounds. We 
perform LSDA+L/ calculations to get a deeper insight. 

Within LSDA-I-C/, depending on the choice of the in- 
put density matrix, we were able to correlate (fill the 
spin-up band and unfill the spin-down band) both band 
complexes in the vicinity of the Fermi level. Contrary to 
other typical cuprates, the ground state is obtained when 
correlating the 'broad' Cu d^z'^^f2 orbital (i.e. (i~2_j^2, 

d\,2_y2, d\^2_f2 ■ filled and d^^2_f.2 ■ empty) with an en- 
ergy gain of about 110 meV (^ 1280 K)'23'] per Cu atom 
as compared to the 'narrow' c?j2_^2 orbital (i.e. d^^^2_f2, 

dii2^f2, dl2_y2 ■■ filled and dl2_y2 ■ empty). This result 
is in accordance with the second scenario described above 
and also with the ID TBM result. Superexchange via 
the apical Cu-0 orbitals is quite surprising and unique, 
though not categorically excluded. Since the apical Cu-0 
bondlength is slightly shorter than the in-plaquette one. 
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TABLE I: Hopping parameters and the corresponding ex- 
change constants from a two-site one-band TBM. The hop- 
ping paths used are: ti [000]-»[100], t2 [00Q]^[110], ta 
[000]-^[110], t4 [000]^[001], ts [000]^[101], te [000]^[iii], 
t7 [000]^[|i|]. The hopping paths are indicated in the inset 
of Fig. [21 The fit has been done for both the in-plaquette 
d^2_y2 and the out-of-plaquette d3j2_i;2 bands separately. 
The parameters are vanishingly small when no value is pro- 
vided. 

the covalency of Cu with the apical oxygens are slightly 
larger than that with the planar oxygens. We can also no- 
tice from the LDA band structure (Fig. [5]) that the d^^-.p 
band is quite narrow compared to d3z^_f2 band. Keep- 
ing this in mind, we can interpret the unique orbitally 
ordered ground state to be the result of the relative gain 
in kinetic energy when correlating the broad Cu d^^2_f2 
orbital. 

Our TBM for the Cu d3z2_f2 orbitals results in very 
short-ranged interactions (Table. H]). Therefore an exten- 
sion of the supercells beyond the second neighbors is not 
necessary. LSDA-I-J7 calculations of differently ordered 
spin configurations (FM and AFM) were performed to 
obtain an effective exchange constant, Jj" by mapping 
the energies obtained from the calculations to a Heisen- 
berg model (neglecting the other J's since they are small 
according to the TBM), H = J2{i j) Jij^i ■ Sj which leads 
to Ef,n - Eaf,n = 2 J2"" | S'l^ , = i. This modcl leads 
to J2" ^ 140 K, in good comparison to the value of 
about lOOK obtained in susceptibility experiments. [iGj 
The slight overestimation of the leading exchange is 
rather typical for this type of calculation and also 
somewhat dependent on the choice of U which is not 
exactly known. We obtain a charge transfer gap of 2.2 
eV which is consistent with the green color of the sample. 



The presence of quite regular octahedra in CuSb206 in- 
troduces an orbital degree of freedom among the eg or- 
bitals as compared to other standard cuprates. Elec- 
tronic structure calculations give convincing evidence for 
the strong competition between the dx2-y2 and d^z^^f^ 
orbitals for the ground state in CuSb206. Correlations 
drive a unique d3£2_^2 orbital ordering pattern in the 
ground state (Fig. [3]). Our TBM and total energy cal- 
culations (mapped onto an effective Heisenberg hamil- 
tonian) are in favor of a ID (d^z^-r^) instead of 3D 
{dx2_y2) scenario following a particular orbital order- 
ing. Our ordering pattern is consistent with the magnetic 
unit cell proposed from the results of neutron diffraction 



3 




FIG. 3: (Color Online) Middle panel: The competition 
of in-plaquette and out-of-plaquette orbitals in LDA for 
CuSb2 06. Both orbitals are partially occupied and comprise 
the Fermi surface. Right panel: Standard picture of high Tc 
cuprates with the inclusion of strong Coulomb correlations U. 
The plaquette orbitals remain half-filled and thereby would 
lead to 3D magnetic order. Left panel: Unique orbital or- 
dering in CuSb2 06 leading to a strong ID magnetic order. 



experiments. [13, This in turn elucidates that the low- 
lying excitations for cuprates are not restricted a priori 
to the standard orbital picture. 

The analysis above has only considered the isotropic 
superexchange interaction. Allowing for spin-orbit cou- 
pling introduces anisotropic terms in the superexchange 
interaction. Many model calculations have been per- 
formed for various cuprates to determine the magnetic 
anisotropy energy, 2J| which scales inversely with the 
difference between the d^2_y2 and 6-3^2 _j,2 energy lev- 
els. Contrary to other cuprates, this energy difference 
in CuSb2 06 is extremely small due to the presence of 
a rather regular octahedron. Therefore, we can qual- 
itatively infer that CuSb206 exhibits a large magnetic 
anisotropy. This anisotropy reduces the phase space of 
the isotropic Heisenberg model and should lead to a sig- 



nificant increase of the magnetic ordering temperature. 
Indeed, if we estimate the Neel temperature (Tjv) as- 
suming an isotropic S = ^ quasi-one-dimensional chain 
within the bosonization method I25ll26l. 



T. 



N 



0.33/cJ'z, Jin 



5.8J 



N 



2 Tn 



we find T^'*^" « 1.5 K which strongly underestimates the 



observed T^''* 



.5 K. Due to the mean field- like treat- 
ment of the inter-chain coupling, T^*^" is expected to be 
an upper estimate. Q The strong underestimation of the 
ordering temperature using the isotropic model implies 
the presence of considerable anisotropy in the system, 
which is naturally explained by our orbital scenario. 

In conclusion, electronic structure and TBM calcula- 
tions in CuSb206 reveal ID behavior from orbital or- 
dering due to a competition between the Cu (ij2_y2 and 
d3j2„f;2 orbitals; strong Coulomb repulsion selects the 
latter (Fig. [3]) . The main reason for such a ground state 
is the underlying local environment. Our results are 
in agreement with the susceptibility and recent neutron 
diffraction measurements for the spin structure. (l6j Per- 
forming polarization dependent x-ray absorption spec- 
troscopy (XAS) should elucidate the choice of orbital 
ordering in CuSb2 06. A strong anisotropy in the ab- 
sorption edge will be observed when using the incoming 
radiation E || a-b plane as compared to E || c. We are 
currently analyzing the results from this experiment. Ef- 
fects of pressure on this system are quite interesting, as 
they can flip the balance between the two competing or- 
bitals thereby restoring the standard plaquette picture as 
in other cuprates. 
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